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Abstract 



CN ' The element of met ast ability is incorporated in the fluctuating nonlinear 

hydrodynamic description of the mode coupling theory (MCT) of the liquid- 
' glass transition. This is achieved through the introduction of the defect den- 

o ■ 

' sity variable n into the set of slow variables with the mass density p and the 



momentum density g. As a first approximation, we consider the case where 



' motions associated with n are much slower than those associated with p. Self- 

consistently, assuming one is near a critical surface in the MCT sense, we find 



that the observed slowing down of the dynamics corresponds to a certain limit 
of a very shallow metastable well and a weak coupling between p and n. The 



^ ' metastability parameters as well as the exponents describing the observed se- 



quence of time relaxations are given as smooth functions of the temperature 

without any evidence for a special temperature. We then investigate the case 

where the defect dynamics is included. We find that the slowing down of 

the dynamics corresponds to the system arranging itself such that the kinetic 

coefficient governing the diffusion of the defects approaches from above a 

small temperature-dependent value 7^. 
64.70.Pf, 64.60.My 66.30.Lw 
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I. INTRODUCTION 



While there has been considerable recent progress in the development of the theory of 
relaxation near the liquid-glass transition, these theories have ignored one of the fundamental 
defining qualities of the problem: These systems are metastable. In this paper, we introduce 
the element of metastability into a framework that connects with mode coupling theory 
(MCT) ||T|. We find a theoretical picture which is closer to the observed experimental picture 
than conventional MCT, since it is freed from the idea of a sharp transition temperature Tq. 

The dynamics of supercooled liquids is characterized by very long relaxation times and 
a rapid increase of the viscosity with decreasing temperature. Thermodynamic quantities 
such as specific heat, compressibility, etc. show a change when the temperature reaches the 
so-called calorimetric glass transition temperature Tg which is defined somewhat arbitrarilty 
by the temperature where the viscosity reaches ~ 10^^ poise. Recent interests, however, have 
focused on the study of the relaxation of supercooled liquids at temperatures well above Tg. 
The reason for this focus has been mainly due to the discovery of the existence of a very 
elaborate sequence of time relaxations spread over many decades of time. This sequence has 
been observed in a series of experiments including light scattering P,^, neutron scattering 
^ and dielectric measurements As depicted for the density auto-correlation function 

in Fig. 1, this sequence of time relaxations is a three-step process with characterizing 
exponents, a, b and (3. After a fast microscopic process, decays as a power law, f + Ait~°- 
followed by the so-called von-Schweidler relaxation, / — A2t^. For longer times t > Ta, 
decays with a stretched exponential, exp{ — (t/rQ,)^}, (3 < 1. 

The mode coupling theory has been very successful in explaining this sequence of time 
relaxations. As originally formulated by Leutheusser 0, MCT is based on the nonlinear 
feedback mechanism where the renormalized viscosity is expressed in terms of the integral 
of the quadratic term in density auto-correlation functions. If we ignore the coupling of the 
density to the energy density, the Fourier-Laplace transform of (j){t) can be written in the 
generalized-hydrodynamics form 
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where Qq is a microscopic 'phonon' frequency and T{z) is the renormahzed viscosity and the 
wavenumber dependence has been suppressed. In MCT it is assumed that T{z) can be 
expressed in the form 

r{z) = do + Ql dt e"' Cn0"(t), (1.2) 

•^0 n=l 

where do is the bare viscosity. The original theory due to Leutheusser corresponds to all 
c„ = except C2. This theory does not lead to the von-Schweidler regime and does not give 
stretched dynamics. Instead, one finds the exponential decay (/5 = 1) following the power 
law decay regime. The importance of the term linear in (ci) in T{z) was emphasized by 
Gotze [|10| as a sufficient condition for obtaining the stretched exponential relaxation and 
von-Schweidler regime. 

It was first noticed in Ref. that MCT can be understood in terms of the fiuctuating 
nonlinear hydrodynamics (FNH). This analysis was then formulated in detail by Das and 
Mazenko |]TT[. One can then use a well-developed field-theoretic methods to determine 

T{z) perturbatively for almost any choice of the free energy. In this formalism, it can be 
understood []T1|, for example, how the linear term in T{z) is generated. Gotze had originally 
introduced the ci term on phenomenological grounds |lO(] . 

Despite its success, the current status of MCT is not without questions and controversies. 
First of all, the scaling picture arising from the experiments of Dixon et al. shows universal 
features beyond those suggested by MCT. In their dielectric measurements, they were able 
to find a scaling curve that fits all the data obtained over 13 decades of frequency range 
for various glass forming liquids and for different temperatures. This beautiful result has 



not yet found an explanation within MCT. As indicated by Kim and Mazenko |T5|, the 
data of Ref. [^] also implies a universal relationship between the von-Schweidler exponent b 
and the stretching exponent (3: (1 + + /3) =universal constant. Although this result 
can be made compatible with MCT , it requires the introduction of additional arbitrary 
constraints. 



Another important question associated with MCT concerns the temperature dependence 
of the theory. Conventional MCT [|T],|^|10|,|T6| views the glass transition as a sharp ergodic- 
nonergodic transition as the temperature T approaches the ideal glass transition temperature 
To well above Tg. According to conventional MCT, the relaxation sequence occurs only near 
T Tq and in particular the von-Schweidler relaxation and the stretching are confined to the 
region T > Tq. This sharp temperature dependence is clearly in contrast with the universal 
behavior of Dixon et al. whose data cover any reasonable choice of Tq. Also, as shown in 



other recent experiments the exisistence of a well-defined temperature Tq well above Tg 
is very doubtful. In this formulation of MCT, the exponents governing the time relaxations 
are temperature independent just as those involved in a second order phase transition. This 
is again in disagreement with recent experiments. In particular, Dixon et al. finds f3 weakly 
temperature dependent, which means, by virtue of the universal relationship mentioned 



above, that b also depends on temperature. Since MCT predicts |]T6| that the exponents 



a and b are related through a temperature-independent relation, this also implies that a is 
temperature dependent. 

Das and Mazenko |TT| , in their study of FNH of dense fiuids, discovered that there exists 
an important nonhydro dynamic correction that cuts off the sharp nature of the ideal glass 
transition. This cutoff corresponds to an extra contribution 7(2;) in Eq. (|1 . 1| ) : 

z + iT{z) 

Although it is now widely believed that the sharp transition is indeed smeared by this cutoff 
mechanism, the question of the exisistence of Tq still remains. There exist recent efforts 
T8| to reconcile the discrepancies between experiments and MCT, which includes the cutoff 



effect into conventional MCT but still assumes a rounded transition around T ~ Tq. This 
attempt, however, ends up with adjusting as many as seven independent parameters to fit 
experimental data. A simpler interpretation of the situation is that MCT, despite all the 
successes, needs to be reformulated at a more fundamental level such that it is not tied to 
the notion of a sharp transition temperature. 



In Ref. [|^, we proposed a model where defects and metastabilty play an important 
role in the glass transition in an attempt to reformulate MCT such that it is compatible 
with experiments. Although this new theory does not provide any clear explanation for 
the scaling result of Ref. 0, it does result in a smooth temperature variation without any 
indication of a special temperature Tq. The various exponents, according to the model, are 
weak functions of temperature. In this paper, our main focus is to elaborate on the detailed 
construction of the model, which was omitted in Ref. and to further elaborate on the 



dynamics of the defects which we treated using a simplifying assumption in the previous 
work. 



The model introduced in Ref. |]T9| is based on the introduction of the defect variable n 
into the set of slow variables consisting of the mass density p, and the momentum density 
g, in the FNH description of MCT ||rT|. The introduction of the defect variable would be 
required in a rigorous hydrodynamic description of crystalline solids, along with the Nambu- 
Goldstone modes associated with broken translational symmetry [^. Although there is no 
broken continuous symmetry involved in the glass transition, we consider the situation where 
this variable plays a role in the transition . The best analogy to the present case would 



be the situation where one considers an order parameter in the disordered state. We do not 
need a microscopic definition for the defect variable here. The only information we need 
here is that the defect density has the usual Poisson bracket relations of a scalar variable 
with the momentum density g, that they are metastable and that they interact weakly with 
the mass density. 

One of the key assumptions in the model is that motions associated with the defect 
density n have a very long time scale compared to that of density fluctuation. This is 
realized in the model via an explicit double well potential h{n) for the defect variable with 
the metastable defect density n associated with the minimum of the higher well. A very 
small diffusion coefficient F^, for n results from a rolling around in a shallow metastable well. 
We find that the coupling between the mass density and the defect variables slows down the 
defect motion further. The coupling also enables the slow dynamics of defects to infiuence 



the density dynamics. Over the significant time period where the defect is trapped within 
the metastable well, the defect auto-correlation function ilj{t) can be regarded as a constant 
while the density auto-correlation function (f){t) displays the relaxation sequence. This is 
essentially the basic mechanism discussed in Ref. that generates the linear term in the 



mode coupling integral. Since one variable (in this case the defect variable) is extremely 
slow compared to the other, we can replace a term like ip{t)(f){t) in T{z) by a constant times 
(pit). Under this assumption, we find that the observed stretched dynamics corresponds to 
a certain self-consistent limit of weak coupling and low activation barrier for the defect. As 
will be discussed below, this limit corresponds to the situation where the coupling energy 
is weak enough not to destroy the metastable defects but still strong enough for the slow 
dynamics of defects to result in the slowing down of the mass density variable. In this 
limit, the parameters describing the double- well potential and the coupling, as well as the 
exponents of the relaxation sequence are self-consistently determined as smooth functions 
of temperature. 

Since the defect variable is diffusive, the defect auto correlation function il){t) decays, 
to linear order, as ~ exp(— 7„t), where in terms of the bare diffusion coefficient for 
n, 7t, = h"{n)V^q^ at the wavenumber q. Self-consistently we find, using the previous 
assumption that il){t) is regarded as a constant, which corresponds to taking the limit 7^, — 0, 
that h"{n) must be small, i.e. the metastable wells become very shallow and broad. As one 
moves into the later stage of relaxation, however, we expect that nonlinear corrections to 
7^ become important. Therefore, we need to consider the dynamics of the defects and their 
coupling to the density fluctuations by including the renormalization of 7^ in the evolution 
equation for 'ipit). In this case, T{z) in Eq. (|1.1|) is not represented simply by the equation 
for (f){t) alone as in Eq. (|1.2| ), since the Ci term is now written as il){t)(j){t) with a nonconstant 
il){t). Instead, we obtain a set of coupled equations for (j){t) and il){t) governing the dynamics 
of both density and defect variables. 

We find, through mainly numerical investigations, that the extended model including 
the defect dynamics, provides considerable self-consistent information on the nature of the 
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defects near the glass transition. In particular, we find that the bare value of 7„ is signif- 
icantly restricted when the observed stretching occurs. In fact the long time scale of the 
defect variable corresponds to the case where the system arranges itself such that 7^ is close 
to some small temperature-dependent value 7^. The system seems to pick out its own defect 
potential and diffusion coefficient. We find that, when 7^, ~ 7^, the time scale associated 
with the defects become much longer than that of the density variable, and thus the basic 
picture obtained using the previous simplifying assumption still holds: The parameters de- 
scribing the metastable wells and the coupling arrange themselves such that the coupling 
enables the slow dynamics of defects to effect the density dynamics without destroying the 
metastable wells. 

We note that the approach we take in this paper is not to explain how the defects bring 
the system to the observed slowing down starting from the microscopic description of the 
defects. Instead, we assume that the system arranges itself to be on the critical surface 
associated with MCT. We then investigate the conditions that the defect degrees of freedom 
must satisfy to be on this critical surface. 

The cutoff mechanism of Das and Mazenko |]TT| will eventually influence the very long 
time dynamics by generating an exponential decay. We note that there exists a mechanism 
|p^ , p!5|| that drives 7(2;) to a small value so that the observed slowing down is retained. 
Although, in principle we can include this effect into the coupled equations, we consider 
here the situation where 7(z) = 0. Another simplification made in this analysis is to neglect 
the wavenumber dependences of correlation functions. In terms of FNH, one can construct a 
wavenumber dependent model as shown in Ref . . Since this model is difficult to analyze 
numerically, we focus here on the wavenumber independent model. 

In Sec. we present a detailed formulation of FNH of simple fluids including the defect 
variable. We then use the well-known field-theoretic techniques to calculate the relevent 
nonlinear contributions to the glass transition. In Sec. |ITI|, the case where the defect auto- 
correlation is a constant is considered in detail as a first approximation. This will reproduce 



the results of Ref. [|I9[. The analysis of the full model which includes the defect dynamics 



is presented in Sec. |^ In Sec. 0, we discuss the cutoff mechanism and the temperature 
dependence of the viscosity in this formulation. In Appendix, we consider the case of a 
general potential h{n) without assuming the particular double-well form. 

II. FNH WITH THE DEFECT DENSITY VARIABLE 

In this section, we formulate the FNH of compressible fluids in detail including the defect 
density n as the additional slow variable. Using well-developed field theoretic methods, we 
calculate nonlinear corrections to the density and defect autocorrelation functions system- 
atically. 



A. Generalized Langevin Equations 

Our starting point is the generalized Langevin equation for the set of slow variables, 
ipa = p(x), 5'i(x), ?2(x), where p is the mass density and g is the momentum density. Here a 
labels the type of the field, the position x, and the vector label i. Following the standard 
procedure described by Ma and Mazenko |23|, we have the equation of motion given by 



'-^-VM-Er.,g^e„. (2.1) 

where F is the effective Hamiltonian and is the streaming velocity governing the reversible 
dynamics and is given by 

SF 

KM = E{^-Mt-, (2-2) 

and {ipaiipii} is the Poisson bracket among the slow variables. The dissipative matrix Fq,^ 
and the Gaussian noise Gq, satisfy 

{Qa{t)Qp{t')) = 2kBTT^p6{t - t'). (2.3) 

The effective Hamiltonian for ip^ is given by 

F = FK + Fu[5p\+F,[5p,nl (2.4) 



where Fk is the kinetic energy: 



Fu[5p\ is the potential energy for the density fluctuation, 5p = p — po, where po is the average 
density, and Fv[6p, n] governs the defect density and its coupling to p. In general, Fu[Sp] can 
be any local functional of 6p and the spatial derivatives of Sp. In particular, one can study 
the wave-number dependence of the structure factor by including the spatial derivatives . 



Recently it has been claimed that by including the wave number dependences one can 
effectively generate the von-Schweidler and the stretching. Analytical treatment in this case, 
however, is very difficult. In the present case, we consider the simple quadratic form which 
corresponds to a wave-number independent structure factor, 

F46p] = J rf3x|(5p(x))', (2.6) 

where A is the fiat inverse susceptibility and wavenumbers are restricted to values less than 
a cutoff A. Fy[6p,n\ in Eq. ( |2.4| ) is assumed to be of the form: 

F4Sp,n\ = J £x [B5p{x)n{x) + h{n{x))], (2.7) 

where we introduced the simple coupling term through the coupling constant B. As discussed 
in Sec. ||, we construct the potential energy h{n) to be a double well potential [^] with the 
metastable defect density n associated with the higher well. We parametrize h{n) such that 
h{n) has three extrema at n = 0, (1 — a)n and n: 

h'{n) = en{n-{l-a)n){n-n), (2.8) 

or upon integrating, 

Hn) = en'[l& " ^(2 - + ^(1 - -)(-)^]. (2.9) 

The parameter e gives the correct dimension for h{n) and describes the overall scale of the 
potential energy. For small positive cr, Eq. (p.9|) represents a double well potential with the 
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global minimum at n = and the metastable minimum at n. We note that for small negative 
0", h{n) in the above representation is also a double well potential. In fact, h{n) is invariant 
under a ^ a' = — o"/(l — a) and n ^ n' = (1 — a)n. Thus, in the following analysis, only 
the absolute value of a will play a role. We also note that if a = 0, then h{n) develops an 
inflexion point aX n = n. It is useful to define two dimensionless parameters x and y which 
characterize the scale of the coupling energy and the potential energy respectively: 



The potential and coupling energy are then completely described by x, y, a and n. We 
consider the coupling energy term in Eq. (|2.7| ) as a small distortion to the shape of h{n). In 
particular, n is shifted to n*[(5/9] determined by 







Sn 

If we expand in powers of 5p, 



= B6p + h'{n*[5p]). (2.11) 

n=n* [5p] 



n*[Sp] =n\l + ai(^) + a^i^-^f + ■■■], (2.12) 
^ Po Po ^ 



and we can easily calculate the coefficients Oi, 02, • • • using Eq. (|2.11| ): 



ai = --, a2 = -(^)(-)^ (2.13) 
ay a ay 

We will later consider the fluctuation of n around this shifted n*[5p\: 5n = n — n*[5p]. 
Eq. ( |2.12|) indicates that, for small |cr|, the coupling energy x must be sufficiently small such 



that 

|ai|<l, lasi < 1 (2.14) 

in order to have a sensible expansion. 

Turning to the reversible streaming velocity terms in the Langevin equations, we assume 



that the Poisson brackets involving p and g are evaluated in the usual way as in Ref. pT 
The new variable n is a scalar quantity as p, thus we assume that the Poisson brackets for 
n have the same structure as those for p. The only nonvanishing elements involving n are 
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Hx),^7^(x')} = -V^[5(x-x'V(x) 
Mx),n(x')} = V^,[5(x-x'Mx)] 



(2.15) 



Since there are no Poisson brackets relating p to n and Tpp = 0, the Langevin equation 
for p is simply the continuity equation, 

dp 



dt 



-Vg 



(2.16) 



The streaming velocity for then follows from Eq. 



n = /<'V[{«-(x).p(y)}^ -H tf(x),.(y)}^ + i:{9-(x).«'(y)}^ 



■ 5F ■ 6F 

-p(x)V^^^-n(x)V: 



■-5pi^) 



"5n(x) 



P 



(2.17) 



Using the explicit form of the effective Hamiltonian, Eqs. (^3 -2.7), Vg can be expressed as 
a functional of and n. In terms of the fluctuation, 5n = n — n*[5p], we have, using 

Eq. (|12D, 



^ = -EV,(ai; + a- 



(2.18) 



where 



a-^5,,{cl{5p) + -A{5pf} 



(2.19) 



is the usual stress tensor with the definitions, 



X 



Co = ^Po(l - — ) 

I.41-^-^il±^ 

^ cry a'^y^ 



(2.20) 
(2.21) 



The bare sound speed Cq and the inverse susceptibility A are modified due to the coupling. 
But the terms depending on the coupling constant are even smaller than |ai| and \a2\ by a 
factor of |a;|. Thus we may neglect them and regain the usual relations: 



cl = Apo, A = A. 



(2.22) 
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The new n-dependent terms in Eq. ( p.l8| ) are given by 



(2.23) 



For our purpose, we need only to keep upto quadratic terms in the fluctuating quantities. 
The Langevin equation for g"^ can then be written as 



dt 



-Ev,K. + ar,)-Ev, 



,9' 



(2.24) 



, 3 P , P 

where Lijix) = — ?7o(|ViVj + 5jjV^) — CoViVj with the bare shear and bulk viscosities, rjQ 
and (o, respectively. For later use, we define the bare longitudinal viscosity, Fq = Co + l^o- 



In Eq. (|2.24| ), the Gaussian noise satisfies 



(ei(x, t)e,(x', t')) = 2kBTLij{^)6{^ - x')5(t - t'). 



The streaming velocity for n is given by 



K = / rfVE{-(x),,^(x')}^ = -E VL(Mx 



(2.25) 



(2.26) 



6g^{^') V ■"^■^""P(x) 
Assuming n is a diffusive variable, we must choose F„„(x — x') = — F„Vx5(x — x') and the 
associated noise S satisfies 



< H(x, t)S(x', t') >= -2kBTT,Vl5{^ - ^)5{t - t'). 



We have the Langevin equation for n. 



dn 
'di 



-Ev.i 



9' 



n- 



+ F,V 



(2.27) 



(2.28) 



p ' 6n 

As in the case of g\ we expand n around the metastable state n*[Sp] and write the equation 
in terms of the fluctuation 6n. Using Eq. ( p.l2| ) and keeping up to quadratic terms in 
fluctuating variables, we have 



|(H = -.(l + -^)V.(I)-V.[(Hfl 

+.f„V^[.y(^) - 2.(i±^)(^)(^) + ,(1 + .)(^)^ 

L n a Pq n n 



(2.29) 



+ 2, 



with the rescaled diffusive coefficient 



J- l : 



V — _o J- V 



(2.30) 
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B. Field-theoretic Formulation 



The Langevin equations, Eqs. ( |2.1(j| ), ( |2.24| ) and ( ^.291 ) can be put into a field theoret- 
ical form following standard Martin-Siggia-Rose (MSR) procedures P^JTB|] . It essentially 
amounts to introducing the hatted variable tpa for each field ipa to enforce the equations 
of motion and integrating out the Gaussian noises to get the quadratic form in ipa- We 



introduce the local velocity field V, where g = pV, to eliminate the 1/p nonlinearities ||Tl 
in Eqs. (|2.24|) and ( p.29|) . This relation is then enforced by inserting the identity 



1 = J VV 6{^-V) = \\p\\ J VY 6{g - pV) 
= IIpII j VYVV exp(^z J dl V(l) • (g(l) - p(l)V(l))) (2.31) 



into the functional integral form, where ||p|| is the functional Jacobian between two delta 
functionals, and dl = d^:s.idti, g(l) = g(xi,ti), and so on. As shown in Ref. the 
Jacobian ||p|| has no dynamical effect, so that it can be neglected throughout the analysis. 
The generating functional is then given in terms of a functional integral over "^a, ^o, where 
= {Sp, gi,Sn,Vi}. Without source terms, it is given by 

Z = J V^V^ exp{-S[^,^]), (2.32) 

where the action 5'[\E', \E'] is given by 

S[^,^]=Jdl ^^,/3"^Li,(l)^j -5^/?-^r^V^5^ (2.33) 

dp 
■dt 



+ z5p[^ + V-g] +tJ2V,[g,-pV^ 



dt 

■I 

+ i5n 



'^^"^ + nil + — ) V ■ V + V ■ [{5n)\] 



dt oy' 

- nT,V^[ay{—) - 2x{ )(^)(— ) + (1 + <y)y{—f] 

n a po n n 

where f3~^ = ksT. As a result of Eq. ( ^^.31] ), we have a polynomial action in \E'q and \E'q,, 
to which we can apply the perturbation theory expansion in a standard way. The nonlinear 

13 



corrections coming from the cubic terms p8| in the action generate the one-loop self-energies 
Ti^p which modify the zeroth order inverse propagators [G^^]"^ through Dyson's equation: 



1 f/^O 1-1 



(2.34) 



C. Linearized Theory 

From the quadratic terms in the action, Eq. (|2.33|) , we can easily read off the elements 
of the zeroth order inverse propagator, [G^^]^^. By inverting this matrix, we get various 
correlation and linear response functions. For later use, we list here a few of these functions. 
In this paper, we will be concerned only with the longitudinal parts of those functions 
and use hereafter the notation p, n instead of (5p, 5n respectively for convenience. In the 
Fourier-transformed space, we have 



G°^(q,u;) = -G^*(q,cj 

1 



Wo' 



ay 



(2.35) 
(2.36) 



where 



Do(q, cj) = po(t^ - g Cq) + iujq Fo 



Wo(q, a;) = (a; + iayTyq^)Do - 



(x + ay f 
ay 



-ujq 



(2.37) 
(2.38) 



The density-density and defect-defect correlation functions are given respectively by 



_^2l3-^q' r.x + ay 



o\ 

1^4 



ay 



) TqUJ q + -r^\Do\ 



Fo|u; + iayT^q^l'^ + Apq{x + ayfV^q^ 



(2.39) 
(2.40) 



Other correlation functions involving n are given by 

2l3~^q^ .x + ay 



Gpni.^..^) = Pon- 



|lVo|2 ay 



UJ To + ayT^po{uj - q Cq) 



(2.41) 



and 
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G° (q,u;) = PoG°..(q,c^) = (-)G° (q,a;). 



(2.42) 



We note that the zeroth order response and correlation functions satisfy the following 
fluctuation-dissipation theorem: 



ay po A 
2/?-i 



ImGVq,u;), 



A PP 



(2.43) 



where a = {p,gi, Vi,n}. 

From Eqs. ( p.35[ - p3^ ), we find that if we can neglect the factor 

{x + ayY 



ay 



ay{l - ax, 



< 1, 



(2.44) 



then we have 



G'°-(q,cu) 



1 



(2.45) 



and 



(2.46) 



L>o(q,^^) 

where 7„ = ayV^q^. The assumption in Eq. ( p.44| ) is equivalent to saying that we are 
considering the small \a\ limit, which will be justified self-consistently later when we consider 
the conditions for the slowing down. Eq. ( p.46| ) represents a standard form for the density 
auto-correlation function. The defect auto-correlation function in Eq. ( p^.45D has a very slow 
diffusive mode, u = —i'jv, since it is assumed that 7„ -C 1. This signals a separation of time 
scales between the density and the defect variables as noted in Sec. |. 

We can achieve a further simplification by looking at the following situation. As one 
approaches the glass transition, we expect that, since the viscosity is getting extremely 
large, one eventually reaches a point where 



Po 



(2.47) 
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Then, from Eqs. (|2.38 



, |2.44|) , the denominator Wq in Eqs. ( p.35|) and ( 



2.36| ) reduces to 



where 7^ = Apq/Tq. It follows then from Eqs. ( p^.39[ - P^¥TD that 

1 T) 

GL(q,^)^-(-) 



213-^ In 



(q,^) ^ - (— — ) — ^ (- .., , + , ..J 



~ — (x + ory) G^„(q, uj) + — ( 



n 



(2.48) 



(2.49) 
(2.50) 

(2.51) 



Po cry 

Therefore, in this approximation, we have the density- defect correlation function as a linear 
combination of the density and the defect auto-correlation functions. 



D. Nonlinear Corrections 



Nonlinear corrections to the zeroth order response and correlation functions in the pre- 



vious section are represented in terms of the self-energies through Eq. ( 2.34 ). For example, 
Eqs. (|2.36| , p.3^ ) are modified to yield very complicated expressions: 

D(q,^) 



W(q,o;) 
1 



W(q,a;) 



[uj + i7^(q, uj)] [p(q, uj)uj + ir(q, uj)q^ 
qn{ci, uj) [qu'iq, uj) + iuj^^Jq, uj)] 



where 



W{q,uj) = [uj + i^',{q,uj)]D{(i,uj) 

n'{q, uj) [uj + igS^^(q, uj)] + p(q, a;)S^^(q, uj) 



-q'^u'{q,uj) 



-^Q^vn^^'^) n(q,a;)[^ -q c (q,^)] -ir{q,uj)q (q,a;)] 



(2.52) 



(2.53) 



(2.54) 
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and 

D{ci,uj) = p(q,^)[a;^ - g^c^(q,^)] + ir(q, + igS^^(q, ^)], (2.55) 

with the renormahzation of parameters given by 

p{q,uj) = po~iJ:^^{ci,uj) (2.56a) 

gc^(q, uj) = qcl + %p(q, uj) (2.56b) 

g2r(q, oo) = q'To + ^S|^(q, uj) (2.56c) 

7^,(q, uj) =lv + «S^„(q, ^) (2.56d) 

qn'{q, uj) = g(l - ai)n + S^y (q, u;) (2.56e) 

gu'(q, ^) = g(5po + (Jen^) + %„(q, ^) 

= (a; + ay)M + s^^(q,^). (2.56f) 

In the above equation L denotes the longitudinal part of self-energy. For example, ^y_y. (q, uj) 
is decomposed into the longitudinal and transverse parts: 



As shown by Das and Mazenko |TT]], the self-energy S^^(q, cu) is responsible for the cutoff 
mechanism. In addition to that, we have self-energies S^^(q, cu) and S^^(q, tu) in this case. 
As will be discussed in Sec. |V|, all three self-energies are related to the cutoff mechanism. 
Although the final stage of relaxation is governed by this cutoff effect, we neglect it in the 
present analysis. 

The separation of time scales observed in the zeroth order correlation functions remains 
valid at high orders if the renormalized version of Eq. (|2.44|) holds, i.e., 



u'{q,uj)n'{q,uj) 



p(q,a;)c2(q,a;) 
In this case, Eqs. ( p.52| , |2.53D give us 



< 1. (2.58) 



n („ A piq,uj)uj + iT{q,uj)q^ 

'^'^ p(q, UJ) [uj'^ — g^c^(q, uj)\ + tL (q, ujjujq'^ 
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Eqs. ( |2.59|) and ( p.60|) are the fundamental equations for the dynamics of density and defect 
fluctuations without the cutoff effect. In the following, we calculate explicitly nonlinear 
contributions to the renormalized viscosity r(q, uj) and the renormalized diffusion coefficient 
7(,(q, u;), which will be expressed back in terms of (^^^(q, tu) and (^^^(q, cj). 



E. One-Loop Evaluation of r(q, oj) and 7(,(q, cj) 

The density feedback mechanism of MCT is realized by calculating the nonlinear correc- 
tions to the bare viscosity in the density correlation function of the form, Eq. ( p.60| ). In the 



present case, however, we have two coupled equations, ( |2.59D and ( p.60| ). Thus we must find 
both kinetic coefficients, r(q, tu) and 7(,(q, cu) to complete the equations. In this section, we 
calculate one-loop nonlinear contributions to those quantities. 



p/^ „ „ P/^ ^p 



(a) (b) 



(c) (d) 



(e) (f) 

One-loop diagrams contributing to 5]^(q, u>) 

Fig. 2 

As noted in Ref. the renormalized longitudinal viscosity in Eq. (|2.56c| ) can also be 
represented in the hydrodynamic limit as 

r(q,a;) = ro-^Si^(q,a;). (2.61) 

The relevant one-loop diagrams contributing to this self-energy are listed in Fig. 2. As 
discussed in Ref. the diagrams coming from the convective vertex {gpVV vertex) in 
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the action Eq. ( |2.33| ) just renormalizes the bare viscosity and are irrelevant to the density 
feedback mechanism. 

We note that the approximation given by Eq. ( |2.51| ) generates two kinds of terms that 
are proportional to GnnGpp and GppGpp, respectively, in the expression for r(q, cj). In fact, 
one can easily see that the GnnGpp type terms come from diagrams (b),(c),(d),(e) and (f), 
and the GppGpp type terms from diagrams (a),(b),(c) and (d) in Fig. 2. Thus we have 



roo r 



(27r)3 

l/(i)(q,k)G„„(k,t) + V^^\ci,k)Gppik,t) 



(2.62) 



Gpp(q- k,t) 



where V^^^ and l/^^) 

are appropriate vertices to be evaluated. In principle, the wave number 



dependence of the vertices can be considered for example, by using the spatial deriva- 
tives of the density fluctuations in the effective Hamiltonian, Eq. (|2.6|) . This involves, how- 
ever, very complicated wavenumber integrals. In this analysis, we consider the wavenumber 
independent case. A closely related approximation to this is that the correlation functions 
can be factorized into wavenumber and time dependent parts [^. We assume here that 
Gnni^l, t) = T{q)ilj{t), Gpp(q, t) = S{q)(f){t), T{q) and S{q) are the fiat structure factors given 

by 



(2.63) 



for g < A and T{q) = S{q) = for g > A, where A is the large momentum cutoff. Integrating 
over wavenumbers for the one-loop diagrams in Fig. 2 with the help of A, we obtain 



T{uj) = Fo + keT 



Qtt^ Jo 



dt e*"*[c/iV^(t)0(t) + ^20 (t)], 



(2.64) 



where the coefficients di and d2 are functions of x, y and a. We find, through detailed 
calculation of the diagrams that 



di{x,y,a) 



-4x(i^)2(2 + — ) 
cr ay 



(2.65a) 



y{2 + 3a) + 2x^ 



.l + a. 



a 



ay 
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d2{x,y,a) = l-4x(^^)(l + — ) 

a ay 



(2.65b) 



+2 



a 



ay 



For later use, we define a dimensionless parameter containing an explicit factor of the tem- 
perature, 



(2.66) 



ApQ Qtt'^ 

The renormalization of 7^ is given by Eq. ( |2.56dD . As in Eq. ( p^.61D , it is equivalent in 
the hydrodynamic limit to 



(2.67) 



Among the nonlinear terms in Eq. (|2.33|) that contribute, at the one-loop order, to the 
self-energy S^(q, cu), those of npn and nnn types contain an explicit factor of which is 
assumed to be very small. The other one-loop diagrams that contribute to E^(q, cu), and 
do not have the explicit F^ factor are listed in Fig. 3. 



(a) 



(b) 



One-loop diagrams contributing to $]:j^(q, w) 
Fig. 3 

Evaluating the diagram (a) using Eqs. ( p. 42]) and ( p.51|) , we find that this contribution 
is directly proportional to 

[x + ayf 



ay 



<1, 



(2.68) 



which is our self-consistent assumption, Eq. ( |2.44|) . The evaluation of the second diagram 
(b) yields 



(2.69) 
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where we have used Eq. ( p.42| ). The second term It in Eq. (|2.69| ) is the integral containing 
the transverse part of the velocity auto-correlation function, G'yy(q, cj). We note that in 
the framework of the flat structure factor, Eq. ( |2.63[ ), there is no coupling between the 



density feedback mechanism and the transverse viscosity. Thus, we assume that this term 
is insensitive to any dramatic slowing down. We obtain flnally that 

n r°o 



7;(q,a;) = 7.-(^)'e dt e'^' m<P{t)- (2-70) 



Eqs. (|2.59|) and ( p.60|) together with Eqs. ( p.64|) and ( p. 701 ) complete the specification of our 
model. 

III. THE MODEL WITH A CONSTANT DEFECT AUTO-CORRELATION 

FUNCTION 

Metastability plays an important role in this model, since defects spend considerable 
amount of time trapped within the well. During this period, we may assume that the defect 
auto-correlation function can be regarded as a constant: 

n 1 

Gr,n{ci,u) = 2ni-f — ^Siu;), or ^(t) = 1. (3.1) 
Po cry A 

This approxiamtion corresponds to the limit where one takes the bare value of 7^, to zero. 
In this case, we note that the nonlinear correction to 7i, also vanishes in the long time limit, 
uj ^ 0, since the integrand in Eq. ( p.70|) is a total derivative when ip{t) = 1. Thus the 



self-consistency of the approximation is maintained. 

The coupled equations ( p.59| ) and ( p.60|) then reduce to a form dependent only on (pit) 



and we can make contact with the standard treatment of MCT [|l|. Indeed Eqs. ( p.60|) and 
( p.64| ) with (|2.65|) are just N = 2 case of the standard MCT equations, ( |1 . 1| ) and ( |1.2| ), where 



only the linear (ci) and quadratic (02) terms are considered. We can identify Qq = q^c^, 
do = g^Eo/p, and 

ci = ^di{x,y,a), C2 = id2{x,y,a). (3.2) 
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It is important to note that the MCT coefficients c\ and C2 are expressed in terms of the tem- 
perature and the parameters describing the metastable potential and the couphng between 
the density and the defect variables. We now assume that the system organizes itself to be 
on the critical surface of MCT, which is described below. This will give relations among the 
metastability parameters. 

The density feedback mechanism associated with the representation of the type Eq. ( |1 . 1|) 
with N = 2 was ffist studied by Gotze |jlO|. According to Ref. |jlO|, there exists a critical line 



in the (ci,C2) space separating ergodic and nonergodic regions, where the nonergodic phase 
is characterized by the existence of the limit / = limj^oo <P{t)i / > 0. To study the critical 
condition in terms of the parameters x, a and ^, we follow a more general discussion on 



such model given by Kim and Mazenko [15]. According to Ref. |15], the glass transition can 



be described by the following three parameters: 



ao = (1 - f)VU) 


(3.3) 


a, = {l- ffV'if) 


(3.4) 


X = \{l-ffH"{f), 


(3.5) 



where V{f) = H(f) — //(I — /) and H(f) = J2iLi ^i/* for a general model containing higher 
order terms in the mode coupling integral. In Eq. (|3.5| ), the parameter A is directly related 
to the exponents of the sequence of time relaxations ||10||: 



r(i-2a) r(i + 26)- ^^-^^ 

The ideal glass transition is approached when both o"o and ci are getting small. Let us 
consider the situation where ui = and the transition is approached by taking ao 0. 
Solving Eqs. ( p.3| - p.5|) for ci and C2, we have 



In terms of the parameters, x, y, a and ^, the critical condition is given by setting (Jq = 0: 
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, / N 2A - 1 , , 

di{x, y, a) = , (3.8a) 

d2ix,y,a) = -^, i<A<l. (3.8b) 

We note that we are mainly concerned with the region where |x| <^ 1 in order to be consistent 
with the fact that we considered the couphng term as a small perturbation shifting the 
metastable state as in Eqs. (|2.12|) and ( p. 141) . We also require the solution be consistent 
with the assumption, Eq. (p.44|) . We find that nontrivial solutions to Eq. ( ^Tsj) satisfying 
Eqs. (|2.14|) and (|2.44|) exist only if we take a — > simultaneously with a; ^ 0, while holding 



^ - (3.9) 



for some constant C. Let us briefly discuss the physical meaning of Eq. ( p.9|) . The small \a\ 
limit is a sufficient condition for a separation of the time scales between the defect and the 
density variables as seen from Eq. ( p.44| ). This indicates that the metastable wells become 



very shallow. If the coupling energy represented by |x| is stronger than the above limit, 
i.e. \x\ ~ (9(1 cr|), then the metastable state given by Eq. ( p.l2| ) does not exist, since the 
condition Eq. ( p. 141) is violated. Thus, a strong coupling energy destroys the metastability 



of defects. On the other hand, if \x\ ~ (9(|cr|'^), we have from Eq. ( p. 651) , 

di = -8y, d2 = l + 4:y. (3.10) 

Since y is always positive, di < 0, therefore the system never reaches the critical surface 
given by Eq. ( p.8|) . Thus, if the coupling is weaker than the limit given by Eq. ( p.9|) , the 
slow dynamics of defects can not affect the density dynamics so that the density variable 
does not slow down. Eq. ( |3.9| ) gives the correct relation between the coupling energy and 
the barrier size of the defects in the case of observed slowing down. In this limit, Eq. ( |2.65| ) 
reduce to 

dJx, y, a) = -8(y + C) - 4( — + 3y + 8C)a + 0(a^), (3.11a) 

y 

d2{x, y,a) = l + Ay + 2{3y + 2C)a + 0{a^). (3.11b) 
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Thus, from Eqs. ( p.8|) and ( |3.11| ), we have the critical condition given in terms of the metasta- 
bihty parameters by 



^ 1, 2A + 1, 11 , , , 

c = y = -,i^.-i)- (3.12) 



Furthermore, from Eqs. ( |3.7[ ) and we can easily see that the limit a — can be 

identified with ctq if 

3{— + 3y + 8C) = 2{l-\){3y + 2C). (3.13) 

Therefore, as indicated in Sec. |, we have the situation where the stretched dynamics is 
associated with the weak coupling {x 0) and low barrier {a 0) limit, which is consistent 
with the separation of the time scales between the density and the defect variables as can be 
seen from Eq. ( p.44| ). The condition Eq. ( p.l3|) , when we use Eq. (|3.12|) , can be interpreted 



as a relation between the exponent parameter A and the temperature represented by the 
parameter i^: 



Xn 2[7 + 2A + V4A2 + 22A + 91]J ^ ^ 



Thus the exponents a and b are given as smooth functions of temperature. (See Fig. 4.) The 
temperature dependence of the parameters y and C which describe the potential h{n) and 
its coupling to p is given by Eq. (|3.12| ). All the temperature dependences certainly do not 



show any indication of a special temperature. In the conventional MCT [0,|16|, the control 
parameter ao is assumed to have a temperature dependence as ctq ~ Tq — T, which is the 
origin of the sharp temperature dependences. In this case, however, is proportional to a 
with weakly temperature-dependent coefficient /(O- 

(To = /(Oa = -^[32/ + 2C]^X{1 - Xy a. (3.15) 

We note that f{S,) vanishes at the lower and upper bounds of ^ which correspond to A = 1 
and ^, respectively. Therefore, the above analysis is not applicable to the region near the 
two end points of ^. Our basic picture of the observed slowing down of the dynamics is 
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that the parameters describing the metastable wells and the coupling, while displaying a 
smooth temperature dependence, achieve the critical limit given by Eq. ( |3.9| ). In this model, 
the transition is actually controlled by the parameter a, which is the size of barrier in the 
metastable potential. 




The exponents a and b as functions of temperature represented by 

Fig. 4 

IV. THE MODEL INCLUDING THE DYNAMICS OF DEFECTS 



The approximation, Eq. ( |3.1| ), made in the previous section enabled us to treat the model 
analytically and to obtain the critical condition between the shallow metastable wells and 
the weak coupling. However, when nonlinear corrections to the bare diffusion coefficient 
come into play, the constraint Eq. ( |3.f| ) must be relaxed and the dynamics of defects should 
be taken into account. In this section, we analyze the full coupled equations, ( p.59| ) and 
( p.60| ) with (|2.64|) and (|2.7CI| ). Since the nonlinear couplings in these equations are highly 



nontrivial, an analytical treatment is very difficult. In the time domain, these equations are 
coupled integro-differential equations: 
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4>{t) + do^it) + ^l(p{t) + t ds Hit - s)0(s) = (4.1) 

Jo 

ij{t) + ^.,^{t) + [ ds G{t - s)ij{s) = 0, (4.2) 
Jo 

where 

H{t) = cMt)<P{t) + C20'(t), G{t) = -U{t)m, (4.3) 

with the initial conditions (j){0) = ip{0) = 1, 0(0) = 0, and the definition C, = (a)^^- ■'■^ 
following, we integrate them numerically. We fix do = = 1 and study the relaxation of 
the system at fixed wavenumber q/A = 0.1 as 7^, and the MCT coefficients Ci, C2 change. We 
use the parametrization Eq. (|3.71 ) for ci, C2 and the value of ^ determined from Eq. (|3.14|) . 

We find that the system depends in a crucial way on the value of 7^,. When the value 
of 7„ is large enough, we expect that the defects quickly diffuse away and only the density 
fluctuations are important. It is clearly seen in the numerical integration of Eqs. ([4.1|) , ( [4.2|) 
for flxed Ci, C2 and ^ in Fig. 5. The defect auto-correlation function ip(t) decays exponentially 
for large 7^,. Thus, in the large 7^, limit, the model essentially reduces to the one originally 
considered by Leutheusser where only the C2 term is present. For smaller values of 7„, 
we expect that the system goes into the regime where the result of the previous section 
applies such that the slow dynamics of defects and the coupling between defects and density 
fluctuations play an important role in stretched dynamics. We can see from Figs. 5 and 
6 that, as the value of 7^ is decreased, ipit) 'Pif) are more and more stretched. We 
can not, however, take 7^, to zero, since for small enough 7^, such that 7^, < 7^ for some 
7^, tpit) starts to increase with time and consequently the model becomes unphysical. This 
surprising result restricts the value of the parameter 7^, to 7^, > 7^. It is clear that il){t) 
and (pit) are stretched most when 7^ is very close to but still larger than 7^. In fact, the 
time scale for ipit) becomes extremely large as 7^ 7^ so that it seems to approach some 
plateau value as t — 00. We will discuss later that only when 7^^ ~ 7^ is the time scale 
of %l){t) much greater than that of 0(t), il){t) is more stretched than 0(t) and the picture is 
consistent with the basic picture obtained in the previous section. 
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In adddition to 7i,, the system also depends on the MCT coefficients ci, C2 such that as 
one approaches the critical surface in the (ci, C2) space, the slowing down of (f){t) is obtained. 
In general, we expect that due to the defect dynamics, the explicit form of the critical surface 
will be different from the one given by Eq. ( |3.7| ) with ctq = 0. In Fig. 7, the relaxation of 
(pit) is shown when we change Ci and C2 but adjust 7^, such that 7^, ~ 7^ for each case. We 
note that in the previous case where ip{t) = 1, the critical slowing down of is achieved 
when ci(A) = (2A — 1)/A^, C2(A) = 1/A^. As can be seen from Fig. 7, for the case of A = 0.6, 
the dynamics of defects influence the system such that at (ci(0. 6), 02(0.6)) = (0.556,2.778), 
(f){t) still relaxes with finite time scales. The critical slowing down similar to the one with 
ipit) = 1 and (ci,C2) = (0.556,2.778) occurs at (01,02) ^ (0.590,2.714) with larger ci and 
smaller 02- We find the similar shift in the critical surface for A = 0.7 case as in Fig. 8. The 
relaxation at (ci (0.7), 02(0. 7)) = (0.816,2.041) for the ip(t) = 1 case is similar to the one at 
(ci,c2) = (0.872,1.902). 

The sequence of the time relaxations of (f){t) can be most easily seen in a \og{—t4>{t)) 
versus log(t) plot. In such plot, the power-law and the von-Schweidler relaxations are com- 
bined to form a minimum with asymptotic slopes of —a and b. The later stage of relaxation 
with a peak is well fit by the stretched exponential (Fig. 9). As 7^, 7^, we find that the 
later stage of the relaxation of ip{t) is also well fit by a stretched exponential. We now use 
two stretched exponential forms 

0(t) = / exp(-(t/r)^), m = 9 exp{-{t/rY) (4.4) 

in Eqs. ( [4.1|) and ( [4.2|) . Since, for 7^ near 7^, the time scale r' for ilj{t) is much greater than 
r for (pit), we can write 

exp(-(t/r')^')^l-(^)^'(^)^'. (4.5) 
We find the following qualitative relations among these parameters from Eqs. ( [4.1|) and ( [4.2[) : 
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Numerically, we find that the second term in the right hand side of Eq. ( |4.7| ) is much smaller 
than the first term. Thus, we can rewrite Eq. ( |4.7| ) as 

1 



(7. - ity 



(4.8) 



which means that r' diverges as {■y^ — 7^)~^ when 7^ 7^. By obtaining the parameters /3, 
P', T and t' from explicit curve fits of and ipit) using Eq. ( [4.4| ) for various values of 7^ 
near 7^, we find some evidence for this relation (Fig. 10). Without the second term in the 
right hand side of Eq. ( |4.6| ), it gives the expression for (3 of the conventional MCT with ci 
and C2 {g = 1): P = — ln(2)/ln((l — Ci)/c2/). However, in the presence of defect dynamics, 
Eq. ( [4. 61 ) represents a more complicated expression for f3 which depends on f3' as well as 
Iv ~ 7^- We find again evidence for this relation from our numerical data (Fig. 11). 



2,0 2.5 



The critical value 7^ of 7^ as a function of the temperature represented by ^ (□). The dashed line is a linear fit. 

Fig. 12 

It is clear from these results that the separation of time scales between the defect and the 
density fluctuations occur only when 7^ is very close to 7^. We find that, as a function of ^, 
the value of 7^ shows a smooth linear temperature dependence (Fig. 12). Under the previous 
assumption that the slowing down corresponds to the a ^ (or 7^ — > 0) limit, we were able 
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to determine the temperature dependence of the metastabihty parameters except for that 
of a itself. This was done by assuming that the system arranges itself to be on the critical 
surface. The present analysis of the model including the defect dynamics determines the 
temperature dependence of the parameter 7^, = ayV^q^ when the system is on the critical 
surface. 



V. DISCUSSION 



It was discovered by Das and Mazenko |TT|] that the nonhydro dynamic correction due 
to Sy^(q, cu) included in the representation like Eq. ( p.GOj ) cuts off the sharp nature of the 
ergodic-nonergodic transition. In the presence of the new variable n, however, we have a 
very complicated expression, Eq. ( |2.53| ) for the density response function. In order to see 
the cutoff effect, let the viscosity V become arbitrarily large. Eq. ( p.53|) then reduces to 

Gpp(q,o.)^ , ^ , , ,. ""^^^ , ,^ , , . (5.1) 



uj{uj + igS^ (q, uj)) + iq S^^„(q, uj) ^^^(q. 



If the self-energies and E^;^^, as well as are set to zero, G^-j;(t) reaches a finite value 
as time t — > 00. The presence of these self-energies makes Gpp{t) decay slowly and thus 
provides the cutoff effect. The evaluation of these self-energies in principle can be done to 
yield complicated expressions in terms of and ipi^)- The analysis including this effect 
thus can be performed but will be a very difficult task even numerically. 

Throughtout this paper, we have concentrated on the time relaxation behavior of the 
system near the glass transition. At this stage, it is important to note the implication 
of our model on the temperature dependence of the viscosity: r]{T). Traditionally, many 
different expressions have been used to fit the experimentally observed ri{T). These include 
the Arrhenius form ~ exp{A/T), the Vogel-Fulcher form ~ exp(i?/(T — Tvf)) and the 
power-law ~ |T — Tq\^^ . These forms are able to fit the experimental data only over limited 
temperature ranges. In MCT, the temperature dependence of the viscosity is given by 

il{T)r^T^{T)^\a^l\ (5.2) 
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where 7 = ^ + ^- According to conventional MCT, (Tq is assumed to be proportional to 
To — T and 7 is temperature independent. Thus, the conventional MCT predicts that the 
viscosity shows a power-law divergence as the temperature approaches Tq. In our model, 
the situation is quite different. The temperature dependence of ctq is given by Eq. ( |3.15| ) 
without any special temperature. More importantly, 7 = 1/ (2a) + 1/ (26) in this model is a 
function of temperature which is represented by the parameter ^. Thus we have 

v{o = \mMor^^\ (5.3) 

where /(^) = |[3?/ + 2C]^A(1 — A)^. Therefore, the temperature dependence of the vis- 
cosity is governed mainly by the behavior of the exponent 7(^) in Eq. ( ^.3|) as a function 
of ^, i.e. logri{^) ~ 7(0- shown in Fig. 13, 7(,^) increases as the temperature de- 
creases. We do not expect, however, that Eq. ( ^.3|) can be used directly to fit the exper- 
imental data, since a and b as functions of ^ are mo del- dependent. Especially the rapid 
increase of log?7(^) near the lower bound of ^ should not be understood as a physical re- 
sult as noted in Sec. |T|. But this analysis clearly indicates that our model is consistent 
with the generic feature of the observed 1]{T) which increases with decreasing temperature. 
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The parameter 7 = 



l/(2a) + 1/(26) as a function of C 
Fig. 13 
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APPENDIX: THE MODEL WITH A GENERAL POTENTIAL h{n) 



In the present work, we have used an exphcit double-well parametrization given by 
Eq. ( p.9|) for the potential energy h{n). In this appendix, we show that the similar analysis 
can be carried out for a general potential h{n) without changing any physical results to the 
lowest order of the perturbation theory. The only information we use for h{n) is that there 
exists a metastable defect density n such that h'{n) = 0. Then the potential energy h{n) is 
described by the following dimensionless parameters: 



/i2 



h"{n)n'' 



At3 



Apl ' - Api 



,2 ' 



(Al) 



and in general, pk = h^^^{n)r)^ j (Apo)? ^ = 2, 3, ■ ■ ■. Thus, the double-well potential given by 
Eq. (p.9|) is a particular case where 



P2 = (ry, /is = 2?/(l + a), P4 = Qy, 



(A2) 



and /ifc = for /c = 5, 6, • • •. 

To the one-loop order, however, only two parameters, /i2 and p^ are relevant, so that the 
coefficients di and d2 in Eq. ( p.65| ) are completely described by x, p2 and ^3. The parameters 
/ifc, /c = 4, 5, ■ ■ ■ generate the terms in the action, Eq. ( |2.33| ) that contribute to the two-loop 
and higher order diagrams. This is basically the reason why the physical situation at the 
one-loop order is well described by the double-well potential, which has two independent 
parameters a and y beside n. 

We note that Eq. (|2.44|) was the self-consistent condition that leads to the separation of 
the time scales between the density and the defect variables. For a general potential it is 
given by 

(/^2 + xf 



IJ'2 



< 1. 



In terms of fi2, and x, the coefficients di and c?2 are given by 



di{x,IJ,2,P3) 



Ai2 



pi 



2(P2 + Ps) + 



pi 



X 



C^2(X,P2,P3) = 1-2^(1 + —) + 
P2 P2 



2(P2+P3) + 



pi 



X 



P2 



(A3) 

(A4) 
(A5) 
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As in the case of the double-well potential, a solution to the critical condition given by 
Eq. (|3.8|) which is consistent with Eq. (|A3|) exists only when we take 



x^O, /i2 0, ^-^C, (A6) 

for some C. This limit corresponds to the case where the coupling is weak and the potential 
gets very fiat near the metastable defect density n. We note that /is must be finite, otherwise 
the model reduces to the Leutheusser model [0], where there is no stretching. This indicates 
that the potential is at least a cubic or higher order polynomial in (n — n) near n in order 
to have the stretched exponential relaxation. In this limit, we have 

di (x, /i2, /U3) = -2/1^(7 - 4/i3 - (4 + S/igC + filC'^)fi2 + 0{iJ,l) (A7) 
d2{x, /i2, = 1 + 2/23 + (2 + 2fi^C)fi2 + 0{fil) (A8) 

As done in Sec. PTl] , we can express the metastability parameters in terms of the temperature 

2A+1 



11 ^22 
/,3 = -(4--l), C=^ (A9) 



£ = - [1 - 3(2A - 1) . 

A2L 2[2A + 1 + V4A2 + 10A + 13]J' ^ ^ 



Because of the particular parametrization of fi2 and fis in terms of a and y, the explicit form 
for the temperature dependence of the parameters for the double-well potential is slightly 
different from that of the general potential. (See Fig. 14.) But, any of the physical pictures 
change from the the double-well potential case. Thus, the double-well parametrization is a 
very good approximation to the lowest order that provides a concrete physical picture that 
the observed slowing down corresponds to the weak coupling between the density and the 
defect variables and the low activation barrier for the metastable defect. 
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The parameter A as a function of ^ for a general (solid) and the double-well (dashed) potentials 

Fig. 14 
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Figure Captions 



Fig. 1. A schematic plot of the sequence of relaxation behaviors predicted by MCT; (a) power- 
law decay relaxation / + Ait~"'; (b) von-Schweidler relaxation / — A2t''; (c) primary relaxation 
and (d) exponential relaxation e 

Fig. 2. One-loop diagrams contributing to S^(q, a;). 

Fig. 3. One-loop diagrams contributing to E^(q,a;). 

Fig. 4. The exponents a and b as functions of temperature represented by ^. 

Fig. 5. The defect auto-correlation function ^p{t) for fixed (ci,C2) = (0.556,2.778) and for various 
7^: (a) 2.0xl0-^ (b) l.SxlQ-^, (c) 1.75xl0-^ (d) 1.725x10-^, (e) 1.7225xl0-^ (f) 1.7215xl0-^ 
(g) 1.721x10-3, (h) 1.72x10-3. We can estimate the value of 7^ as 1.72x10"^ < 7^ < 1.721x10-3. 



Fig. 6. The density auto-correlation function (j){t) for fixed (01,02) = (0.556,2.778) and for various 
7^; (a) 2.0x10-3, (b) 1.8x10-3, (c) 1.75x10-3, (d) 1.725x10-3, (e) 1.7225x10-3, (f) 1.7215x10-3, 
(g) 1.721 x 10-3, and (h) 1.72 x 10-3. 

Fig. 7. The density auto-correlation function (^(t) for various (ci,C2); (a) (0.467,2.944), (b) 
(0.488,2.905), (c) (0.522,2.842), (d) (0.556,2.778), and (e) (0.590,2.714). 7„ is adjusted for each 
case such that 7^ ~ 7^. (f) (p{t) at (ci, C2) = (0.556, 2.778) for the case where ip{t) = 1. 
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Fig. 8. The density auto-correlation function for various (ci,C2); (a) (0.719,2.285), (b) 
(0.779,2.135), (c) (0.816,2.041), and (d) (0.872,1.902). 7^ is adjusted for each case such that 
7^ ~ 7^;. (e) (p{t) at (01,02) = (0.816,2.041) for the case where V(i) = 1- 

Fig. 9. The soUd lines are log(-i^(i)) versus log(t) plots for {01,02) = (0.467,2.944) and for (a) 
7^ = 1.7225 X 10~^ and (b) 7^ = 1.75 x 10~^. The dashed lines are the stretched exponential fits, 
(p{t) = /exp{-(i/r)'^}, where (a) / = 0.44, p = 0.804, r = 5859 and (b) / = 0.44, jS = 0.875, 
T = 4021. 

Fig. 10. The best-fit values of [T'r(l + ^)]-i for various values of 7^ (□). The dashed line is a 
linear fit. 

Fig. 11. 2~^/f^ versus ,f- (77)^ plot (□). The dashed line represents a linear fit. 

Fig. 12. The critical value 7^ of 7t, as a function of the temperature represented by ^ (□). The 
dashed line is a linear fit. 

Fig. 13. The parameter 7 = l/(2a) + 1/(26) as a function of the temperature represented by ^. 

Fig. 14. The parameter A as a function of ^ for a general (solid) and the double-well (dashed) 
potentials. 
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